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INTRODUCTION 


In this report mathematical physics arguments useful for lightning 
discharge and generation problems are pursued. The idea of using a more 
mathematical approach than heretofore is based on the prospects for 
clarification of the physical phenomena involved which in many cases is 
obscured by the multiplicity of factors at work. 

The first section treats a soliton Ansatz for the lightning stroke 
including a charge generation term which is the ultimate source for the 
phenomena. In this way the number of functions required to completely specify 
electric and magnetic fields, charge and current densities is reduced to a 
minimum. For purposes of simplification it is supposed that the ionization 
channel radius is independent of time. 


„ The 3e c on d section establishes dynamical, electrical and thermal equations 
for a partially ionized plasma including the effects of pressure, magnetic 
field, electric field, gravitation, viscosity and temperature. From these 
equations is then derived the Non-Stat ionary Generalized Ohm's Law essential 
for describing field/current density relationships in the ionization channel of 
the lightning stroke. Arguments are then given for the essential participation 
of ionic generation processes in the "exponentially'' increasing current density 
and charge density which develop during the stroke. 

.. t Th u thlrd 3ection deals with, the discharge initiation problem and argues 
that the ionization rate drives both the convective current and electric 
displacement current to increase "exponentially" but that because of relative 
saturation of the former compared to the latter the convective current is 
unable to "relieve" the electric field which eventually increases to breakdown 
unleashing the lightning stroke. In this section the non-linear term of the 
on-Stationary Ohm's Law is retained without approximation so that the temporal 
development of the lightning discharge may be precisely formulated. 


.. fourth 3ec tion deals with the statistical distributions of charge in 

the thundercloud preceding a lightning discharge. Defining centers of positive 
and negative charge, centers of generation and recombination as well as centers 
of current efflux and influx for the cloud it becomes possible to statistically 
characterize the development of the cloud dipole moment and the relative 
velocity between the centers of charge as functions of time. 


The fifth section contains some physical comments on the stability of the 
pre lightning charge distributions and the use of Boltzmann relaxatlonal 
equations to determine them. Also the argument for aircraft providing a 
lowered impedance path for the stroke is given subject to the additional effect 
or field enhancement factors such as aircraft curvatures. 


*NASA-ASEE 1984 Summer Fellow 


tfts-n wz? 



MATHEMATICAL PHYSICS APPROACHES TO 
LIGHTNING DISCHARGE PROBLEMS 


LIST OF SYMBOLS 


q charge per unit volume (net signed charge in a unit of volume) 

q + absolute (non-negative) quantity of positive charge per unit volume 

q_ absolute (non-negative) quantity of negative charge per unit volume 

J net signed charge crossing unit area per unit time: electrical current 

density 

J + absolute (non-negative) quantity of positive charge crossing unit area 

per unit time: electrical current density of (absolute;) positive charge 

J_ absolute (non-negative) quantity of negative charge crossing unit area 

per unit time: electrical current density of. (absolute) negative charge 

note: The directions of the current densities are given by the vectorial 

average of the velocities of the individual charges crossing unit 
area in unit time; with this convention the above definitions of the 
magnitudes of the current densities are completed vectorially 

G net signed charge generated per unit volumd per unit time (includes ion 

production and recombination) 

G + absolute (non-negative) quantity of positive charge generated per unit 
volume per unit time 

G_ absolute (non-negative) quantity of negative charge generated per unit 

volume per unit time, 

c light speed 

t elapsed time 

t distance unit of elapsed time - product of light speed and elapsed time 

V del operator with scalar components representing partial 

differentiations in the coordinate directions 

u magnetic permeability of medium (plasma) 

z electric permittivity of medium (plasma) 

E electric intensity field vector 
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e(r ) 
h(r) 


magnetic intensity field vector 
charge (absolute) of electron 

valence of ion, i.e. non-negative integer multiple of e which is 
absolute ionic charge: for electron or singly- charged + ion Z=1 

diffusion coefficient (one third of square of step length per mean free 
time in a random walk) 

drift velocity in a random walk 

unit vector orthogonal to axis of cylindrical coordinate system 

unit vector in direction of circumferential (e) increase in cylindrical 
coordinate system 

unit vector in axial direction of cylindrical coordinate system 
operator of partial differentiation in direction of u 

r 

operator of partial differentiation in direction of u 

0 

operator of partial differentiation in direction of u 

z 

radial component of magnetic field intensity 

circumferential component of magnetic field intensity 

axial component of magnetic field intensity 

radial component of electric field intensity 

circumferential component of electric field intensity 

axial component of electric field intensity 

fraction of light speed at which soliton (pulse) propagates 

radial dependence of axial component of electric field intensity 

radial dependence of circumferential component of magnetic field 
intensity 


k reciprocal channel size constant (of ionization channel) 

K Q (kr) modified Bessel function of second kind 
I Q (kr) modified Bessel function of first kind 

radial component of current density vector 
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circumferential component of current density vector 

axial component of current density vector 

absolute number of particles per unit volume 

absolute number of positively charged particles per unit volume 

absolute number of negatively charged particles per unit volume 

absolute number of uncharged particles per unit volume 

valence of positively charged particles (averaged for aggregate) 

valence of negatively charged particles (averaged for aggregate) 

mass density: total mass of particles in a unit volume 

mass current density - momentum density: total vectorial momentum of 

all particles in a unit volume 

number of neutral particles generated per unit volume per unit time 
(including recombination and ionization) 

Eulerian derivative operator for differentiation with respect to time 

while following moving portion of plasma D = 9 + v.V 

t t 

magnetic induction field vector 

vectorial average velocity of positively charged particles 

vectorial average velocity of negatively charged particles 

vectorial average velocity of uncharged particles (averaged over small 
portion of plasma) 

mass of positive ions (averaged over small portion of plasma) 

mass of negative ions (averaged over small portion of plasma) 

mass of uncharged particles (averaged over small portion of plasma) 

vectorial momentum transfer to positive ions per unit volume per unit 
time 

vectorial momentum transfer to negative ions per unit volume per unit 
time 

vectorial momentum transfer to neutral particles per unit volume per 
unit time (momentum transfer excludes contributions to pressure and 
viscosity) 
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partial pressure due to positive ions 

partial pressure due to negative ions (including electrons) 

partial pressure due to neutral particles 

viscosity for positive ions 

viscosity for negative ions 

Viscosity for ne.utral particles 

gravitational potential 

arithmetic meaii viscosity 

viscous averagd velocity of small portion of plasma 
resistivity of [plasma 

reciprocal valence weighted mass of charged particles - 1/(z + m +Z m ) 
geometric mean mass of charged particles - (m + m ) 1/2 
valence weighted arithmetic mean mass of charged particles 
average valence of ♦ and - ion 3 
average squared valence of + and - ions 

harmonic mean mass of charged particles weighted for valence and 
concentration 

valence weighted concentration average for charged particles 

geometric mean valence of charged particles 

mass weighted valence difference parameter 

Laplacian operator 

electrical collisional frequency 

vectorial cyclotron frequency for positive ions 

vectorial cyclotron frequency for negative ions 

absolute temperature (Kelvin) 

thermal current density 



k thermal conductivity 

S power per unit volume radiated 

c + specific heat at constant pressure per unit mass for + ions 

c_ specific heat at constant pressure per unit mass for - ions 

c Q specific heat at constant pressure per unit mass for neutral particles 

S + power radiated per unit volume by + ions 

S_ power radiated per unit volume by - ions 

S Q power radiated per unit volume by neutral particles 

T + temperature of + ions. 

T_ temperature of - ions 

T q v temperature of neutral particles 

k + thermal conductivity for + ions 

<_ thermal conductivity for - ions 

•c thermal conductivity for neutral particles 

D + diffusion coefficient for + ions 

D_ diffusion coefficient for - ions 

D q diffusion coefficient for neutral particles 

X mean free path length (t » mean free time) 

v mean free speed 

Y + collisional frequency for + ions 

Y_ collisional frequency for - ions 

Y q collisional frequency 

Cp specific heat per unit mass averaged for +, - and neutrals 
<m> arithmetic mean mass of particles 
<D> average diffusion coefficient 

<Y> average collisional frequency 
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Boltzmann constant 
conductivity 
volume of region 
surface of region 
volume average 
surface average 

unit normal vector for surface 
time dependence of charge density 
time dependence of electrical current density 
time dependence of electric field' 

vectorial electric field intensity generated by + charges 

vectorial electric field intensity generated by - charges 

total absolute positive charge 

total absolute negative charge 

center of positive charge 

center of negative charge 

portion of total volume over which net production of + ions occurs 
portion of total volume over which net recombination of + ions occurs 

center of generation for + charges 
center of recombination for + charges 
center of generation for - charges 
center of recombination for - charges 

portion of closed surface over which current due to + charges is efflux 
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S + portion of closed surface over which current due to + charges Is influx 

S_ + portion of closed surface over which current due to - charges is efflux 

S portion of closed surface over which current due to - charges is influx 

R* + center of efflux of current on surface S due to + charges 

R*_ center of influx of current on surface S due to + charges 

R* + center of efflux of current on surface S due to - charges 

R*_ center of influx of current on surface S due to - charges 

N + total number of ♦ charges generated per unit time in V 

N_ total number of - charges generated per unit time in V 

A total absolute quantity of positive or negative charge (equal) 

v pel relative vectorial velocity of centers of positive and negative charge 

M vectorial dipole moment of thundercloud 

B. Boltzmanriian derivative operator for differentiation with respect to 

. time while following states of small portion of plasma in phase space 

B t “ V + ^ ' V R + * * V 

dts differential volume element in locational space 

rv 

dt~ differential volume element in velocity space 

note: occasionally the same symbols have been used for different variables but 

the difference in meaning should be clear from the difference in context 
and in location on the list of symbols: one should especially note the 

difference between t mean free time and distance unit of elapsed time; k 
reciprocal channel size constant and the Boltzmann constant; viscosity 
and resistivity n 
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1. Soliton Ansatz for the Lightning Stroke with Charge Generation Term 

The basic objective of this section is to derive a mathematical framework 
in terms of which observational and experimental data can be used to identify 
and interpret physical features of the propagation of fields, charges and 
currents during a typical lightning stroke. By characterizing the field and 
current density components in terms of a small number of scalar functions it 
should become easier to relate input data to the actual development of charge, 
currents and fields. 

The starting point is the set of Maxwell equations modified by the 
Inclusion of a charge generation term G describing the number of charges (of 
valence Z) generated per unit volume per unit time. The absolute charge 
density q is taken to be positive corresponding to the absolute value of 
charge per unit volume. Then in principle each of the following equations 
after (1.2) is doubled with a separate version for positive and negative 
charges and for the currents and fields associated with them. The actual 
signed charge density q and current density J are related to these unsiged 
quantities by 


q - Q + - Q_ (1.1) 

J - J + - J_ (1.2) 


but for conciseness the generics q, J will be used for typical q + , or 

q_, J_ corresponding to G f , G_. Generally the inclusion of a G term means 
that charge conservation will be violated for ♦ and - species separately but it 
will be expected that overall conservation holds for the net charge given by 
( 1 . 1 ). 

Turning now to the Maxwell equations one has using t a ct and Gaussian 
units 


V x 

E - -u3 H 

T 

(1.3) 

j - 

c (V x H - c3 t E) 

(1.U) 

q - 

eV • E + ~ f Gdt 

(1.5) 


C J 0 
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H - 0 

(1.6) 


Non-Conservation of charge then follows by adding the partial time derivative 
of (1.5) to the divergence of (l.il) 
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(1.7) 


3 1 ^ + V • J » ZeG 

At the outset it may be noted that the assumption 


J » -DVq + wq 

converts (1*7) to a Diffusion Equation with drift 


( 1 . 8 ) 


8 fc q - DAq t w-Vq + qV*w + ZeG 


(1.9) 


S 1 S e .SI 6 .! f 3 “f 1 ? 1 ? 8 Stroke lfc 13 to be e ^eoted that the DAq terra 
'■ b }p b ! fc n randomizin S element in (1.9) will be small compared with the 
drift term w*Vq producing an often small random deviation in direction of the 
stroke (whose direction is principally determined by w). This point which is 
of passing interest will not be further pursued here. 

in ZurTJilr 6 Equ ? tions (1..5). (1,6) are to be solved 

in cylindrical coordinates with a soliton Ansatz for the case of no 0 

dependence. Components in the 9 direction will, however, be retained The 

thMp S ni°H?M° r CUrl and divergence ln < r » fi.-z) cylindrical coordinates under 
these conditions are recalled to be unaer 


_ 3 (rH ) 

V . H . -u^H, . u 9 (3 z H r -3 r H z ) . u z S—J- 


V • H - 13 (rH,) OH 
r r r z z 


( 1 . 10 ) 

( 1 . 11 ) 


al!ial directions 3 ^ ^ UnU VeCt0r3 in the radial » circumferential and 
By virtue of (1.11), (1.6) may be written 


3 r {rH r> - 3^-^) 


( 1 . 12 ) 


So that the magnetic intensity components 
function V(r, z,t) by 


may be expressed in terms of 


a 


rH 

r 


3 7 


z 


(1.13) 
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(1.14) 


-rH » 3f 
z r 


If V(r,z,T) • rh(r)F(z-ST) one has 


H p - h(r)F'(z - St) 


u . " Crh(r)] , F(z - St) 
z r 


( 1 . 15 ) 

( 1 . 16 ) 


This suggests an Ansatz fop H. of the form 

O 


H 0 « h(r)H(z - St) 


(1.17) 


and a corresponding Ansatz for E of the form 

z 


E - e(r)E'(z - St) 
2 


( 1 . 18 ) 


The cylindrical components of the Faraday Induction Law (1.3) are 


~3 E - -p3 H 
z 0 H t r 


3 E - 3 E - -m3 H 
z r r z t 0 


(1.19) 

( 1 . 20 ) 


-3 r (rE 0 ) - -u3 t H z 


(i.21) 


which upon substitution from (1.15), (1.16), (1.17), ( 1 . 1 3 ) become 


3 z E Q - -uSh(r)F (z - St) 

3 z E p - e(r)E(z - St) - uSh(r)H'(z - St) 


( 1 . 22 ) 

(1.23) 


? 3 r (rE 0 ) ‘ tC!?<r ) 3 , F »(z - St) 


(1.24) 
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Clearly, (1.22) 13 satisfied by 


Eg - -u5h(r)F(z - £x) 


in which case (1.2^) is automatically satisfied. According to (1.23) 


3 E - e(r )E(z - 5 t) ♦ y5h(r)H(z - $ T ) 


which is satisfied by 


E r - e(r)E(z - ^t) + u£h(r)H(z - £x) ( 1t 

Collecting the components of electric and magnetic intensities one then has 


E « e(r)E(z £x) 


E 0 - -uCh(r)F(z - 5 t) 


Ep - e(r )E(z - £x) + u£h(r)H(z 


. [EM£>2' F <z . tT , 


H 0 - h(r)H(z-ex> 


Hp » h(r )F(z -£x) ( 1>3 

For the generation term G in (1.5) and (i.7) one takes 

G(r,z,x) - £g(r) P (z-^x) ( 1 . 3 ; 

where r is an as yet unspecified soliton function (not the Gamma function). 
In terms of r the integral of G can be written 


Gdx - Cf(z) - r(z-£x)]g(r) 


(1.35) 
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It should be noted that (1.5) assumes no generation before time t= 
(1.35) In (1 .5) one has 

1 " 3 z E z * r 3 r (pE r ) + g(r) Cr(z) ' r (z’^T)] 
or substituting from (1.28), (1.30) 

3 - e(r)E , ’(z-CT) ♦ E(z-^t) ♦ u5 H(z-*x) 

p* If 8(r) [r(z) - r(z-5t)] 

At time t-0 this becomes 

f ■ «(de"(z) . tcalicill e(2) , „ 5 [rhicul „ (2) 

For zero initial charge density q-0 

e(p)E* ' (z) * CisjriL E(z) . u5 i£!!iEn'„ (j , . o 

and it must be satisfied with r independent of z. This can be 
accomplished by 

k 2 e(r) - t re ' (r)] '- „ c Cfll(r)l ' 
r r 

and 

e"(z) - k 2 E(z) + k 2 H(z) - 0 

r 

The first equation of ( 1 . ^0) implies 
e (r) + - k 2 e(r) - 0 


0. Using 

(1.36) 

(1.37) 

( 1 . 38 ) 
(1.39) 

(l.iJO) 

(1.H1) 

( 1 . * 12 ) 
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which is satisfied by the modified Bessel functions I (kr) and K (kr) of 
which K (kr) exhibits exponential decay as r •* « wRile I (kr) °becomes 
exponentially infinite as r ♦ «. Thus K (kr) is the betteR solution and one 
may take 


e(r) - A K Q (kr) (1.43) 

The second equation of (1.40) is satisfied by the choice 

I 9 

u£h(r) - e (r) - AkK Q (kr) (1.44) 

while (1.41) is satisfied by 

-k 2 H(z) - e”(z) ♦ k 2 E(z) (1.45) 

Substituting these results into (1.28), (1.29), (1.30), (1.30. (1.32), 

(1 .33) one has 


E z “ AK 0 (kr)E (z ~SO 


E 0 - -AkK^(kr )F* (z-£t) 


E r ' "k <o (l<r)E " (z ' 5T) 


H z - TT k 0 ( 'x-Wz-M 


H 9 ' l^k K o (kr) * e"(2-«t)3 


H r ' ;K <kr)F ’ (z ‘ 5,) 


The components of the Ampere Circuital Law (1.4) are 


J /c - L d (rHj - e9 E 
z r r 0 t z 


J a /c • (J H • 3 H ) ■ ei E 
0 z r r z t 0 


J /c ■ -3 H - c3 E 
r z 0 t r 
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(1.46) 

(1.47) 

(1.48) 

(1.49) 

(1.50) 

(1.51) 

(1.52) 
0.53) 
(1.54) 



so that substitution from (1.28), 


( 1.29), ( 1.30), (1.3D, (1.32), 


(1.33) yields 


V c - Ct) + e5e(r)E(z -£r) 


J fl /c - h(r)F (z - 5 t) + LliliLU' 


F(z - lx) 


- euC 2 h(r)F' '(z - £ T ) 


J r /c - -h(r)H(z - «T)-*.ccCe(r)E(z - 5t) + uth(r)H(z - Ct)] 
for the current densities. 

Substitution from (1.43), 0,44) and (1.45) then yields 
-AK (kr) 


J /c 
z 


[k 2 E(z-£T) + (1-€U^ 2 )e"(z-^t)] 


AKjkr) 

^ - [k & F(z-CT) + (1 -eu^‘ ; )f"(z-5t)] 


J /c ■ - ° *T-r _ \ J. ^2 ,f 

9 


AkK (kr) , 

J r /C * CE (z-5t) + "~o( 1 -eu5 2 )E ' ' ' (z-^t) ] 

k 


while (1.39) applied to (1.37) results in 
p " If C r (z) “ r(z-Si)]g(r) 


which for g(r) - e(r) - AK Q (kr) is 


e " FT K Q (kr)[r(z) - f(z-^x)] 


(1 .55) 


(1.56) 

(1.57) 


(1.58) 

(1.59) 

( 1 . 60 ) 


( 1 . 61 ) 


( 1 . 62 ) 


4 ...T 1 ? 6 3 ° lution mi ght be further simplified by taking F-0 if experimental 

justification excluding circumferential currents could be found. This would 
leave only E(*- 5t ) and r(z-5x> in the solution considerably simplifying the 

stlroke atl R n i- 0f fl ® ld data wlth the Propagation properties of the lightning 8 
stroke. Return strokes are easily handled in the Ansatz by replacing f by -£ 

result^F-E^ls obt^nT?? It be , a Pr ° Per fractlon of u 8ht speed. Another 
vector is set equaJ ?o zer^ component of the Poynting 
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It seems clear physically that a .constant pulse shape can hardly be 
maintained Indefinitely as it propagates along its ionization channel. 

Lt d ° e ! n0t follow that ln fche ^ imes considered the pulse shape will 
change very much and in any case it is primarily important at this stage to 
discern the way in which the various vectorial components depend upon the E 
. «# ^notions. In this way it is hoped that measurements will eventually 

clarify the entire phenomenological sequence of events of the stroke. 


2. Plasma Dynamics and the Non-Stationary Generalized Ohm’s Law 

The purpose of this section is to systematically construct the dynamical 
electricai, and thermal equations describing the behavior of a typicar P ortion 
nL? 1 ? ma c °^ tainin 8 positive and negative ions (including electrons) and 
_ rr r f 1 , p f rticl ® s entirely free of mathemati cal approximations. This h» 

effected by means of suitably defined averages. T he electrica l equations are 

SSdT S T U- relating llZrull 

.I* 1 ® J' ormu J atlon win include effects of partial pressures due to ♦ and 

oni ?« a ? d #.? e ^ Pal partlcle3 * external electromagnetic fields, a gravitational 
potential field and viscosity terms. ■ 

There will naturally still remain physical approximations and statistical 
wnria ions concerning the extent to which such averages can reasonably 
represent the parameters entering into the formulation. 

^ The restriction here to an isotropic model is not serious since the 

iriPfl? U !n 3 <!I a f be reapil y generalized to stresses as in the transition from an 

? d ® inli? 1 ? t0 an an } aotr °P lc viscous fluid. The tensor extensions possible in 
the model are generally easily handled and will be ignored here. 

First relevant averaging parameters are defined and then equations for 
mass, charge, momentum, and thermal transport and generation are established. 

2a. Definitions and Formulations * 

The notation adopted associates a + subscript with positive ion 3 , a - 
subscript with negative ions (including electrons) and a zero subscript with 
neutrals. The ions are supposed to be of valences Z + and Z . 

We define 


i$ns l<jns rjieutrals 
n - n + + n_ + n Q - number density 

q - (n + Z f - n_Z_)e - charge density 


( 2 . 1 ) 

( 2 . 2 ) 
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mass density 


(2.3) 


P “ nj. + n m + n m - 
+ + — . oo 

J - n m nmv+nmv-pv (? h\ 

mass current density (defines v) 

(also equal to momentum density) 

J ■ Z + en + v + - z _en_v_ - electrical current density r (2.5) 

so that 

V ' m + 3 t n + + m - 3 t n - + m o 3 t n o (2.6) 

7 • J ■ m 7*(n v.) + m 7*(n v ) + m 7*(n v ) (2 7) 

m ▼ ▼ — —— OOO • t / 

3 t q “ Z + e3 t n + " Z - e3 t n - (2.8) 

V • J - Z + eV*(n + v + ) - Z_eV*(n_v_) (2.9) 

We assume no net mass generation and no net charge generation 

3 t p + 7 * 5 m “ 0 “ 3 t q + 7 ' 5 (2.10) 

No net mass No net charge 

generation generation 

With G + , G_, G the generation rate for + ions, - ions and neutrals (the 
latter is initially negative and finally positive in a transient ionization- 
recombination sequence) one has 

3 t n + + 7 * “ G + generation rate of + ions (2.11) 

3 t n _ * 7 * ( n _0 “ G_ generation rate of - ions (incl. electrons) (2.12) 

3 t n o + 7 * (Vo 1 " G o g enerati ®h rate of neutral's ( 2 . 13 ) 
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so that the ionic generation rates are related to the neutral rate as shown 
below (2.16) (2.17) 


ra^C, +mG +mG » 0 no net mass generation 
+ <• - - oo 


(2.14) 


Z + G + - Z_G_ - 0 no net charge generation 


(2.15) 


Consequently, 


G 


+ 


G 


-m Z 

o ~ G 

m_Z f + m + Z_ o 


m_Z + + m + Z 


- number of positive ions generated 
per unit volume 
per unit time 

G q - number of negative ions (or electrons) 
generated per unit volume per unit time 


( 2 . 16 ) 


(2.17) 


The electromagnetic fields are regarded as externally imposed, it being 
supposed that the particle motions are sufficiently random as to produce 
effective cancellation of the fields generated by individual particles except 
when their motions become so highly correlated as to produce well defined 
resultant fields. 

The Eulerian derivative D. ■ 3. + v • V of the mass current density J 

t t m _ 

(mass/area time - momentum/ vol. ) and of the electrical current density J are 

given by (2.18) and (2.22) 


D. J - nrD (n v ) + m D.(n v ) + m D. (n v ) 
tm + + + - t - - otoo 


( 2 . 18 ) 


JxB»J + xB- J_x B 


(2.19) 


J + - en^Z + v + 


( 2 . 20 ) 


J - en Z v 


( 2 . 21 ) 


D fc J - eZ + D t (n + v + ) - eZ_D fc (n_v_) (2.22) 

By solving for v + and v_ in the expressions for the two current densities J 
and J these quantities may be expressed in terms of them and of the velocity 
of neGtral particles v q as shown in (2.26) and (2.28) 
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m 


en + Z + v + - en_Z_v_ - J 




eZ + [n + m + v + + n_m_v_ 


J - n m v 
ra ooo 


eZ 


en (Z m + Z m )v - eZ.J - m,J - en m Z v 
+ + + m + o o + c 


eZ J - m J - en m z v 

+ ro + oo+o 

en_(Z + m_ + Z_ra + ) 


en (Z m + Z m )v - eZ J + m <J - en m Z v 
+ + - + + - m . - oo-o 


eZ J + m J - en m Z v 
- ra oo-o 

■ ^ — — — — — — — — — — — — 

en + (Z + m_ + Z_m + ) 


m n v 


[eZ J - m J - en m Z v ]m 
+ ra + oo+o- 


e(Z + m_ +'Z_m + ) 


m + n + v + " 


[eZ J + m J - en m Z v ]m 
-ra - 0 0-0 + 

e(Z + m_ + Z_m + ) 


The momentum densities for the 3 species satisfy the equations 


> t (n + m + v + ) - en + Z + E + en + Z + ( £ ) x B - Vp + - n + m + V* + T + 


+ n + Av + + - - V(V • v + ) 


D t (n_m_v_) - -en_Z_E - en_Z_(~)x B - Vp_ - n_m_V<|> + T_ 


+ n_Av_ + — V(V • v_) 


D t (n o m o v o ) * n o Av o + T 7(V * v o ) ' ? P 0 ' n o m o 7; t ) + ^o 


( 2 . 23 ) 

( 2 . 24 ) 

( 2 . 25 ) 

( 2 . 26 ) 

( 2 . 27 ) 

( 2 . 28 ) 

( 2 . 29 ) 

( 2 . 30 ) 


( 2 . 31 ) 


( 2 . 32 ) 


( 2 . 33 ) 



The terms T + , T_, T represent momentum transfer rates per unit volume to 
♦ Ions, - ions and neutrals due to collisions. <j> is gravitational potential 
and n . n and n are viscosities. Addition yields 

+ “ O 


D t <J m - D fc (pv) • (J x B)/c - 7p - p7<J> + 3n^A < v > + n A (7 • <v>) (2.3^) 

under the assumption that the net tranfer of momentum for the whole system is 
zero. Thus T + ♦ T_ + T q » 0. <v> is a "viscous average" velocity defined by 

n + v^ + n_v_ + n 0 v 0 " 3n A <v> with n A arithmetic mean viscosity. 

The electrical current densities for the ions satisfy 


D t (en + Z + v + ) 


2 2 
e n + z ; - 

2 _2 
e n * Z + 


eZ . 



(v + X 

B) - — Vp 

■ 11 ■ Cj ▼ 

m + 

m + c 

eZ + 

eZ + 


eZ T 

+ ST n+Av+ 

3m. 

n + 7(7 

• * ST 


en + Z + V<j> 


D fc (en_Z_v_) 


2 2 
-e n Z c 


m 


2 _2 
e n Z 


E - 


eZ ' 


m c 


eZ 


eZ_ 

(v_x_B) - jjj— - Vp_ - en_Z_7<f> 


eZ 


V - — n Av + r— - n 7(7 • v ) + — T 

m ■" * ■< m — m — 


3m 


m 


These may be subtracted to yield 


V 


v m. c ' m m / 


m + m 


ae 

( z i*. 

c 

' m . 

eZ_ 


m_ 

Vp_ - 




m + m 


ae ZZ n m / Z Z \ 

+ -oo/ + - \ - - 

- I — + — j v xB 

c \ m m_/ o 

e *- e S /Z + _ z - - \ 

“ — Vp_ - - — Vp + e(n Z - n Z )7<J> + e ( — T T ) 

m_ m + + - - + + \ m + + m_ -J 


x B 


(2.35) 


( 2 . 36 ) 


♦ eA 


"kiii . Lizil . • „ c, rw* . fiijh 

_ m + m_ _J 3 \ m f m_ / 


(2.37) 
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where a Is defined by a(Z t m_ + Z_m + ) - 1. A resistivity n may be defined 
in terms of the momentum transfer terms for ions thus 



n - resistivity (not to be confused with viscosities n + , n_, n Q ) 

By introducing the following means, the electrical current density equation can 
be made more concise 


m„ - m rn 

u ~ + 


m. 


( Z + m + + Z_m_ v 
Z~ + Z / 


<Z> - ^ Z + + Z_ ^ /2 


2 

<Z > 


( z * * z - ) 


/ 2 


2 2 
n.Z n Z 


1 /«V + v 

— - l — — + — — ) /(n Z? + n Z ) 
m_. \ m. / + + - — 


"h 


n ft - ( n + Z^ ♦ n_zl)/(zf * zf) 


2 < Z > n. 


"h 


2 2 

riX 

m + * 


Z G ■ Z - Z + 


(Z + m_ + Z_m + ) 


(2.39) 

(2.40) 

(2.41 ) 

(2.42) 

(2.43) 

(2.44) 

(2.45) 

(2.46) 

(2.47) 
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(2.48) 


Z^m_ Zjn + 

Z <5 " m + ’ m_ 

L - Laplaclan operator, 7 - del operator 
In these terras, one has 


V 


2 2 

2e < Z > n / 
m H 


e 2 Z*n ra _ / Z_ 2 + \ 

( x B) + e ( __ Vp Vp J ♦ e(n_Z_ - n + Z + 

m G C 

•¥]•!'(••[¥•'•€]) 


(I - nJ) + 


2„2 
e Z G P 
2 

m G c 


cteZ 


(v x B) + 


J x B 




(2.49) 


Defining the conductivity a by 
1 

a - - 

n 

and the electrical collislonal frequency r by 
2 2 

2e^ < 1 > n A 

r - n 

H 

and Introducing the vectorial cyclotron frequencies 


(2.50) 


(2.51) 


Z + eB 


m + c 


(2.52) 


fl 


Z_e§ 

m c 


( 2 . 53 ) 
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The second and third terms may be replaced so that 


V 


r(oi - J) + 


otZ, 


eZ_m f n + _ _ eZ + m _n_ 

■ (Q + x vj — ~ ~ (G_ x v_) 


m 


[r < 5 * 5 -> - < J * 5 .>] * * Cf - r t *.) 


/ z + n + v + Z_n_v_\ 

♦ e(n Z - nZ x )V* + eA l ) 

" * + N m . m/ 


+ r V 


\ nr m / 


(2.54) 


The simplest non-stationary case is with B 
pressure and gravitational gradients. 

Then (2.54) becomes 


0 and negligible viscosity, 


D fc J - r(oE - J) • (2.55) 

as the generalized Ohm's Law which reduces to E - nJ or J » oE for 
- - - z _ z + 

D fc J - 0 - B - V<J> and — Vp _ - — Vp + with n + - n_ - n Q - 0 (2.56) 

(2.55) will be considered in more detail at the end of this section. 

In terms of g the vectorial acceleration due to gravity (2.54) can also be 
written 


a E - 


V 


2 2 

2e* < Z* > n. 


V + 


m H Z G P0 


2 < Z 2 > n A m 2 c 


(v x B) 


™H Z 6° 


2ce<Z >n. 


. Z G n o m o m H° r. 

J X B r- (v x B) 

2 < Z^ > n A mQC 0 


2e 


m H° / z _ z + \ m„a 

2 l nT Vp - " F 7p - ) + 5 (n_Z_ - n + Z + )g 

< Z d > n A ' m - • / 2 e < Z 2 > n. 
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V ■ 

2e < Z 2 > n. 


m H° 

6e < Z 2 > n , 


Z f n + V <- z _n_v_ 


m 


V 7 


z + r ’+ v + z _n_v_ 


m. 


m 


(2.57) 


pressure) p specific heat/mass at constant 

PV t T * V . j„ - E . J/O . „ t y t . * v** • 7,(7 • V 


n v • Av + 


n_(v_ • 7)(7 • v ) 


n (v • V) (V • v ) 
+ n v • Av + -2 — 2 o 


o o o 


+ S 


S is the energy radiated per unit volume per unit time 
temperature grldfenT"" 1 ° Urre3t de " 3Uy 13 take " *> Proportional to the 


(2.58) 


J 0 - - < V T 


(2.59) 


one may write for the S^peciero^pfrticle^ 6 generalized 0hm ' s Law - Hence, 


V" + c + 3 t T v - 7 *(k + 7T + ) 


n_m_c_3 t T_ - 7 • ( k_7T ) 


en + Z + (E«v + )/ c + n + v *Av + ~ (v *7)(7 • v ) + s 

(2.60) 

-en_Z_(E*v_)/c n_v_-Av_ + ~ (v_«7)(7 . v_) + S 

( 2 . 61 ) 
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n oVo 3 t T o " V ,(< o VT o ) ' Vo'^0 + f (v o .V)(V.v o ) + S o 

where T, T_ and T are the temperatures of + ions, - ions and neutral 
particles. 

A plasma temperature may be defined by 


n.m c T +nmcT + n m c T 
♦ + ▼ + — — — — , 


OOOO 


n.m c. +nmc +nmc 
+ + + ___ ooc 


Then a specific heat at constant pressure c y may be defined by 


n.m.c + nmc + nmc - pc 
+ + + - - - ooo p 


where 


n.m ♦ n ra -t-nm - p 
* * - - oo 


Hence 


W. T . ' n -"-°- T - * Vo°o T 0 ■ P V <2.6l 

Diffusion coefficients may be defined by 

n + m +c + D “ < + > n _ m _c_D = k , n m c D - < (? 6i 

ooooo U " U| 

in terms of the thermal conductivities k + , < q for the 3 species. 

s lnce diffusion coefficients in general are obtained from a random walk model 
with mean free time t and mean free path A one has 


D ■ r 


where v 


Hence 


( 2 . 69 ) 
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and since 


mv 2 2 
— - | kT 


( 2 . 70 ) 


one qbtains 


D - !Sll . !<I 
m mY 


( 2 . 71 ) 


with colllaional frequency y. Substi t-nt-i n« 

the expression for pc kT one has 8 3 result for ea °h species Into 


o o -o'V'o pC p kT 


Y + n + m + c + D + + Y_n_mfc_D + y n m 2 c D 

o o o o c 

If <m> denotes the arithmetic mean mass with 
3<m>«m+m + m 

o 

define a statistical average collisional frequency < y > b y 

- y > i&ak i * T - n - m - c - * y - n - m2 ° 

m > (n + m + c + n_m_c + n m c ) 

ooo 

or 


( 2 . 72 ) 


( 2 . 73 ) 


(2.7H) 


V * n * n *°* * T - n -"- 0 -* Y o" 0 "o° 0 ■ < V > < m > 


pc. 


30 that an average diffusion coefficient < D > 


becomes 


< n v I*WA * Y _ n _ m fc_D + y n m 2 c D 
^ D > - O o o n o 

^ Y > < m > pc 

P 

Hence the expression for pCp kT yields 


< D > - 


kT 

< m > < y > 


( 2 . 75 ) 


( 2 . 76 ) 


( 2 . 77 ) 
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Returning now to consideration of the Non-Stationary Ohm's Law (2.55) and 

approximating the Eulerian derivative D J by 3 J one has 

t t ^ ~ 


3 fc J - r(oE - J) 


(2.78) 


Then (1.7) can be partially differentiated with respect to time to yield 


3 t q + 


V*O t J) - Ze3 t 0 


(2.79) 


into which (2.78) may be substituted to yield 
2 - 

3 f q + ToV*E - rv*j - ZeZ.G 

Substitution from (1.5) and (1.7) then yields 

3 t q + r 3 t q + r 2 q “ rZeG + Ze3 t G + f l Gdt 


( 2 . 80 ) 


(2.81) 


in which the charge generation 
Inhomogeneous equation. For 
equation 


terms act as forcing functions in the above 
G-0 one obtains the corresponding homogeneous 


a t q H + r3 t q H + r q H ' 0 
which is satisfied by 
„ 

0„ - q,e 


f or 


(2.82) 


(2.83) 


. rx + — 

€ 


<2.81 0 


which is seen to have no roots with negative real part since 

21 - -r ±"Vr 2 - 


(2.85) 


corresponding to exponential decay or oscillation under 
decaying envelope in (2.83) depending upon whether 


an exponentially 
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(2.87) 



or 

JJo 

r < ^ (2.88) 

In any case the argument indicates that decay of charge is inevitable 
unless G ^ 0. Thus nonvanishing charge generation is a sine qua non for an 
"exponentially" increasing charge density. 

In the next section examples of increasing charge will be derived without 
neglecting the non-linear term in (2.55). 


3. Lightning Stroke Initiation Problem 

The purpose of the section is to attempt to characterize the mathematical 
features of the initiation problem in sufficient detail to clarify the physical 
theory of the phenomenon and provide a framework for the design of experiments 
and observations of the development of lightning discharges. 

The starting point is the 3 et of Maxwell equations modified by inclusion 
of a charge generation term G which is a function of position R and time x 
- ct signifying the number of charges (of valence Z) generated per unit volume 
per unit time. 


V x E - -m3 H 


(3-D 


J - c(V x H - 

e 9 E ) 

(3.2) 

q - eV-E + — 
M c 

f Gdt 

(3.3) 


J o 


V-H - 0 


(3.4) 

Elimination of E 

between the time derivative of 

(3.3) and the divergence 

of (3.2) leads to 



3. q + V-J « ZeG 

(3.5) 
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In this q 1 3 meant to be the absolute quantity of a v, 

to 1 the° 1 h' ne i'1 Ule t U the 00rre3 P°"< 11 "8 current density and 8 G refer 8e only 

vSw n ZT °: ZTZZTl °a h n a d r8 1 TillT 8eneraUd dS 

these (single-signed) charges. Hence^ (L^pen^ Tm^ions 


' •> 

3..q. + 17*3 - ZeG 


. 

t + + + 

(3.6) 

: ' * 

\ t 

9 f q + 7*J - ZeG 


! ; 

t - 

(3.7) 


with the actual (signed) charge density q and actual 
given by 


current density J 


q - q + - q_ 


J - J - j 

■ T • 


(3.8) 

(3.9) 


sped f icat ioidboth 6 redaU1 " 8 “ at “ lth tbl3 

the Oausfddrgenc; 5 Theore m a y?e?5s y ' 00nne0tad V ° 1Ume V b0Unded by * « 


<3 t q> v + V X "•‘j > s - z e <G> 


( 3 . 10 ) 


notes tha£ for V ° 1Ume ^ 3Urface avera S es respectively. One 


<3 t q> v ’ ”f < n-j> 3 


(3.11) 


chAr^in* 13 ^ £ 7 33ible to have an increase in the absolute quantity of 
charge in a region together with an efflux of ciFri^t uni M « . «, JL fMirr 

The normal vector n is taken to poinF^ard (i^rom V) on sT P ° SitiVe * 
prevlIusS5 e der?ved? qUaU ° na °" e al3 ° addS the "“"-Stationary Ohm's Law 


D fc J - T(oE - J) 


(3.12) 
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which may also be written 


q3^J + (J*V)J - qroE - qTJ 

If a solution of this partial differential equation be sought which is 
separable in q, J and E one has 

_ a 

q(R, t) - q(R)a(i) 

J(R, t) - J(R )6( t) 

E(R, t) - E(R )Y( t ) 

( 3 . 13 ) becomes 

A A | A A A A A A 

cqJaB + 0 ( J • V) J =■ qEfoaY - qJfctB 


( 3 . 13 ) 


( 3 . 14 ) 

( 3 . 15 ) 

( 3 . 16 ) 


( 3 . 17 ) 


or 


A A I A A A, A A 

qJ( TB+cB )a + B (J*V)J - qEroaY 

which can be satisfied separately in time and position provided 
(B + ~ B )ct ■ B^ = otY 

AAA 

and (absorbing constant factors in q, J, E) 

A A A A A A 

rqj ♦ ( J*7) J - qETo 
In other words 

Y - B ♦ p g' 

TB 2 

« “ — r 

TB + cB 


( 3 . 18 ) 


( 3 . 19 ) 


( 3 . 20 ) 


( 3 . 21 ) 


( 3 . 22 ) 
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t 

For a temporally Increasing current density 6 >0 and (3.21) indicates 
that the field will increase faster than the current density. Also 


i - 

i + 

I 


i 


i 


» f 0 t t 

V * 6 + 

For 

8 - (e XT - 1) k - W k 
with w'- X(W + 1) 

• it - 1 ' 1/ 1/- 1 

0 - kW K W - Xk(W + W ) 

e" - Ak(kW k-1 + (k-1)W k “ 2 )X(W+1 ) 
or 

i i o I , u - 1 It-? 

0 - X*k[kW K + ( 2k- 1 )W K ' + (k-1 )W K *] 


(3.23) 


(3.2H) 


(3.25) 


(3.26) 


whence 

(3.27) 

(3.28) 

(3.29) 

all of which are initially zero. 

For times x >> 1 /X (t >> 1/Xc) one has the asymptotic formulae for the 
temporal dependence of current density, electric field, charge and displacement 
current 


Y a w k + X£k w k-1 


i 


Y - Xk |(1 ♦ ^)W k t I 1 + 


r + Ac(2k-1) j.jk-1 , Ac(k-I) W k-2~J 


and by (3.22) 


a - W k+1 / (l + W + 


B 


Xkx 

e 


(3.30) 



(3. 3D 
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- ^ kT / l i ^ \ 

a e / ^ 1 - — ; 

Y 1 - Xk (l + ) e XkT 


so that Ak can be chosen sufficiently large that the ratio of displacement to 
convective current is large. 

In the presence of "exponential" forcing from the ionic generation term G 
both convective current and displacement current will increase but eventually 
the convective current will be "relatively saturated" with respect to the 
displacement current which thus provides the physical mechanism for increase of 
the field to breakdown. It should be clearly understood that the displacement 
current referred to here is monotonlcally increasing not oscillatory as in most 
electromagnetic applications. 

After breakdown occurs, it leads to the development of an ionized channel 
with enhanced conductivity so that it then becomes possible for larger 
convective current to flow which ultimately "relieves" the high E field and the 
lightning stroke proper is over. The disturbance in charge balance can, 
however, lead to return strokes and further effects. The constants r and X 
must be determined experimentally. 

Naturally the increase in q, J and E will not occur everywhere but 
only in locations subject to (3*20). Applying the Gauss Divergence Theorem to 
( 3 . 20 ) one has 


f < qj | < (J-n)J > 3 - f < oqE > y ( 3 . 34 ) 


or 


2 A ^ AAA 

v < >3 ■ r < 0 (qE - J) > v (3.35) 

which suggests that the larger the discrepancy from the static Ohm’s Law the 
less likely the discharge region is to be spherical since the sphere has the 
smallest surface S for a given volume. 

It should be particularly noted in the arguments of this section that the 
results (3.24) , (3.29), (3.27) and ( 3 . 28 ) for the current density, charge, 
field_and displacement current take full account exactly of the non-linear 
(J*V)J term in (3.13) the Non-Stationary Ohm's Law appropriate for this case. 
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1 t. Tho Separate Charge Distributions in a Thundercloud 

In this section the statistical properties of the separate ♦ and - charge 
distributions in a thundercloud will be investigated and relations involving 
the cloud dipole moment will be established. 

Writing (1.5) for the separate charge species one has 


z + e f T 

q+ - ~ / G + dT + eV*E 
' o 


(4.1) 


Z_e + t 

q_ » J G_dr - eV*E 


(4.2) 


for the positive and negative absolute charge densities. 

Using the Gauss Divergence Theorem this leads to total charges 



(4.3) 


(4.4) 


In terms of these (positive) total quantities of positive and 
one can define 


negative charge 


Q R V* V 


'4.5) 


as the center of positive charge and 



(4.6) 


as the center of negative charge . 

The rates of change with time of the total charges of the two species is 
found from (4.3) and (4.4) to be F 



d t Q - " z _ e J G _dV - ,/n.^.dS 


( 4 . 3 ) 


^velocities of the + and - charge centers can be found fran ( 4 . 5 ) 


R * " d t R + " Q J R3 fc q + dV - R + d t (J?.n Q + ) 

(4.9) 

R - “ d t R - " Q_J R3 t q_dV - R_d t (£n Q_ ) 

(4.10) 

Substituting from the modified continuity equations 


3 t q + “ Z + eG + " 

(4.11) 

3 j.q_ ■ Z_eG_ + V* J 

(4.12) 

(4.9) and (4. TO) become 


• Z + e r - 

R * " 57 J RG + dV - 5 + J R(V-J + )dV - R + d t (^n Q + ) 

(4.13) 


Vt 



z - e r - ,f. . 

q— J RG_dV ♦ ^ J R( V»J_)dV - R.d (An-Q_) 


(^ . 14) 


30 that if V ++ is the portion of total volume V over which net production 
of + ions is occurring and V + _ is the portion of total volume V over which 
recombination of + ions is occurring one has V - V + v and one mav 

ripflnA J 


■ vf J 

r- 

| R|0jdV 

(4.15) 


V ,. 

5° . J- 1 
*- v 

f R|GjdV 

(4.16) 

+ - ~ 

' V - 


as the center of generation for + charges and the center of recombination for 
♦ charges. Similarly 


*?+ * \T“ J R|GjdV 
- J V - + 

f R|G_|dV 
-- J V 


are the center of generation for - charges and the center of recombination 
for - charges. The terms in (4.13) and (4.14) containing the divergences of 
current densities can be written in terms of surface integrals by the Gausa 
Divergence Theorem. Thus 


J"r ( V • J ) dV - <j) 


R(n»J)dS 


(4.19) 


so if the portion of S in (4.13) corresponding to current efflux is S and 
the portion corresponding to current influx is S + _ one can define * + 




I. 


R | (n*J + ) |dS 


++ 


(4.20) 


35 



(4.21) 


I 


R*_“ s J R I l ds 
+ - S + - 


as the centers of current efflux and current influx due to positive charce. 
Similarly in (4.1 4) 


R* + - 5 1 R | (n • J ) | dS 

- + s - + 

R L- § | R | (n*J_) | dS 


(4.22) 


(4.23) 


In terms of (4.15), (4.16), (4.20), (4.21), (4.22) ,• (4.23) one has from (4,13) 
and (4.14) 


- z + e -r -r ( S ++ RL- S + r[ ) 

K “ q 7 ( v ++ r ? + - VR?_) q- “ R.Mln QJ 


+ t' 


(4.24) 


Z e 


xG 


R_ " q~ (V_ + R_ + ~ V R__) + 


G . S„rL) 


R_d fc ( £n Q_) 


(4.25) 


Returning now to (4.7) and (4.8) one notes that according to (1.4) 

1 


fn-^E dS - s n*VxHdS-^^|"n*JdS 


and 


J n*7 x H 


dS - 0 

over a closed surface so 


(4.26) 


e (p n*3 E dS - - I 


^ n*3 fc E_dS » 


(4.27) 

(4.28) 
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where I + and I_ are the currents flowing put of V due to .+ and 
charges respectively. Thus (4.7) and (4.8) become 


Q + - d Q + - Z + eN + - I + 


(4.29) 


d t Q_ - Z_eN_ + I_ 


(4.30) 


where N + and N_ are the total number of charges + and - generated in V 
per unit time. i j 

In case of net charge zero - Q_ - Q and the relative velocity of the 
+ and - charge centers is given by 


so 


V rel‘ R . - R - 


« rei - w° - w 1 - 


Q(R + - R_ )d t ( in Q) 


(4.31) 


(4.32) 


with 


W G - Z + e(V ++ R G + - V + _R G _) - Z_e (V_ + R_ G - V__R G _) 


W 1 - (S t+ Rj + - S + _R*_) + (S_ + R* + - S R*_ ) 


(4.33) 

(4.34) 


The dipole moment M of the cloud is 


M - Q(R + - R_ ) 


(4.35) 


and its time rate of change is 


M - Q V re i + ” d t Un Q) (4.36) 

The condition for constant separation of the charge centers is V , = 6 hence 
(4.34) implies rel 

M d^ffcn Q) - W G - W* (4.37) 
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Thus the averaging arguments to establish the statistical structure of the 
cloud charge distributions is established. 


5. Physical Arguments Bearing on Lightning Discharge Problems 

According to Earnshaw's Theorem electric charges cannot be maintained in 
equilibrium by electrical forces alone so it is important to recognize that 
static models have a limited use in explaining lightning and the development of 
thunderclouds capable of generating it. Actually one can start with static 
models in unstable equilibrium and augment the electric and gravitational 
forces by aerodynamical flows giving at least temporary stability to the 
structures while the charge concentrations develop in the cloud. The proper 
tool for the temporal development of these + and - charge distributions is 
the Boltzmann equation 


B ,.P - 3..P ♦ vV p + a*V p - S p (5.1) 

1 1 R y 1 

where 


p(R, v, t) dt dT (5 2 ) 

R v 

is the probability that a charge is located within dx R of R and with a 
velocity within dr of v and 5 p is the 30 -called "collisional 
derivative" which is not actually a derivative at all but an integral giving 
the rate of change of p due to collisions. Undqr some conditions which are 
probably adequate for the problem of explaining the physical basis of charge 
separation in the cloud in mathematical detail the 6 p may be replaced by Yp 
to form the Relaxatlonal Boltzmann equation 1 ^ 

B t P “ y P (5.3) 

In the stationary , coll is lonless case (5.3) is easily solved in terms of the 
canomical ensemble for which 


P 


P o e 


-E/kT 


(5.4) 


The problem of charge separation is clearly neither stationary nor 
collisionless so (5.3) must considered as it is and its solution is obviously 
somewhat more complex than (5.4). Since (5.2) must satisfy the normalization 
condition 


If p dx dx 
R v 


(5.5) 


'For negative Y, 5.3 corresponds to decay to a null state; if it is desired to 
decay to a final state p Q , then p should be replaced by p-p . 
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for all time it is clear that Y cannot be taken to be constant since this 
will lead to a probability density with an exponential factor e™ in the 
solution so the (5.5) cannot be satisfied. Thus (5.3) must be studied in 
generality to rigorously establish the properties of the separately signed 
charge distributions and it is intended to pursue this direction in a 
subsequent publication. Here we merely note that it follows from (5.3) that 
the locations of maximum probability density correspond to the positions of 
null resultant force for each charge species and the physical explanation of 
the charge separation process must be found in the fact that the null force 
altitudes are different for positively and negatively charged particles. The 
three types of force term most relevant for the altitudk segregation process 
are the electric field force, the gravitational force and the average force due 
to collisions. The electric field force is altitude dependent and oppositely 
directed for the two charge species. The gravitational force is altitude 
dependent via the altitude variation of the acceleration due to gravity and 
will reflect the mass difference of the two charge species. The coins ional 
drag force will depend on the altitude variation of temperature and 
concentration which will also play a role in the therma . biasing of fracture 
processes in riming. Although the null force condition for the two charge 
species may be fulfilled at separate altitudes it cannot be expected that such 
a separation of charges will remain stable for any appreciable time without the 
influence of a vorticlal flow which would tend to keep the separated charge 
configuration within its low pressure region. 


The average fair weather electric field of the Earth is directed downward 
and decreases monotonlcally .with increasing altitude. Typical values are 120 
volts/meter at the surface, 25 volts/meter at 3 km. and 6 volts/meter at 10 km. 
Below 3 km. the field is sufficiently strong to accelerate electrons upward out 
of the layer of atmosphere closest to the Earth. The acceleration of positive 
ions is less because of the mass difference between electrons and positive ions 
and the net result is to establish an excess of positive charge in this region. 
However the most important features of the situation are not the static ones. 
There is a continuous upward transport of electrons and a continuous downward 
transport of positive ions to the negatively charged Earth under fair weather 
conditions. This is the observed Sky-Earth current. Far above it is to be 
expected that there will be an escape of electrons into the thermosphere 
analogous to the escape of hydrogen since the electric fields are presumably 
negligible and the electron gas is even less dense than hydrogen. This is of 
little relevance to events in the troposphere except as a consistency condition 
for the general upward .flow of electrons. 


With the advent of pre-storm conditions there is a disruption of the 
normal adiabatic lapse rate so that the lower troposphere becomes considerably 
cooler than it was under fair weather conditions. This is due to enhanced 
mixing between layers at different altitudes. In this way the lower 
troposphere attains a temperature intermediate between the fair weather surface 
temperature and the tropopause temperature. This quasi-isothermal layer 
interferes with the vertical transport of heat leading to the development of 
thermal vortices with updrafts near the surface and it also interferes with the 
transport of electrons upward and the transport of positive ions downward so 
that each of these species tends to collect at characteristic altitudes with 
the positive ions usually concentrating at higher altitudes because they are 
subject to much higher drag in the thermal updrafts. There is a very large 
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size difference between the electrons and the positive ions. Again it should 

statiS and^hat^ !? arge 3eparatlon P^ess is dynamical rather than 
fnn ^ ii d ® pends u P° n the aerodynamic flow to remain stable for any 

arromnffi V ength ° f tim ®* After the se P ara tlon of charges becomes y 

*®°°* pli . 3hed °? e can e *P ect that the local effects of the charge maxima will 
ehario i" ® lect ™ n ? being repelled upwards and downwards near the negative 
charge center at altitudes above and below it while positive ions will be 

aid e below r n! Ued UPWaPdS and d0Wnwards near the Positive charge center above 


. nH should be clear that continuous transport of electrons away from Earth 
° " tlnuou3 transport of positive ions to Earth under fair weaker 
nSff!! nS t, Cann0t 5® raaintalned indefinitely if the Earth is to retain its 

Jefult of C the 8 U t rhJ!Iin rder f °H t0 be Preserved one expects that the net 

f lightning ground strokes in thunderstorms will be to return a net 
negative charge to the Earth so that a consistent fair weather electric field 

aeonmnn a h n i a i ne ?‘ WUh regard to individual strokes this could be 
accomplished by transporting negative charge to Earth or positive charge from 

°n de n storm conditions charge separations of 2 to 6 kms. occur with chanze 
°f tbe order of ten3 of coulombs and these generate breakdown ^ 
f*®* d3 3 " Iega '' 0lts/meter in dry air; 1 megavolt/meter in clouds). Since 

breakdown is observed to occur it is a simple matter to compute 


E - k Q 2 /r 2 - 10 volts/meter 


(5.6) 


with k - 9 x 10* 


so 


Q/r ~ .01 


(5.7) 


for TxoTZI bre3kdow " there mu3t be roughly 10 coulombs of charge of each sign 
for each kilometer of separation to generate such a field. 8 

o The !! air W ® ather surface field of about 120 volts/meter corresponds to a 
surface charge density of about a nanocoulomb/meter 2 . This yields an Earth 

^ 0 - C me£?rs? f "SX m ® gacoulomb — the total areao^Jhe Earth of 5 

x 10 meters . With an observed air conductivity of 2 x 1 0~ l Vmpi-pr =. 

Sky-Earth current density of about *.« plcoamperes/meter* Is obutneS tance 
the total current over the whole Earth is 1200 amperes according to the 

fh rm H° f ? hm,S LaW * ThlS mU3t then be compensated during lightning 
strikes in thunderstorms by much larger current acting for much shorter time! 


6. Summary and Conclusions 

The basic purpose of the soliton Ansatz used in the first section «« 

* 9 number ° r run = t ‘°'» in terms of which the 3-oomponent eUMrlc Held 
intensity, convective current density, magnetic field intensity and the scalar 
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electrical charge density and Ionic generation may be described. The intention 
is that the functions left unspecified in that section be determined by 
experimental data. Once they are specified the full functional dependence of 
all components would be determined. Also the relationship of the input data to 
the full characterization of the vectorial fields would be precisely .specif led 
and the directional properties of the lightning strike fields clarified. 

In order to understand the processes underlying the initiation of 
lightning discharges it is essential to have a physically reasonable, 
mathematically manageable formulation of the relationship between the electric 
intensity field and the electrical current density. The purpose of the second 
section is to give a systematic plasma dynamic derivation of the nonlinear, 
nonstationary generalization of Ohm's Law appropriate for the plasma of the 
ionization channel. The simplest case (2.55) suitable for the lightning 
discharge problem involves an electrical collisional frequency in addition to 
the conductivity characteristic of the stationary form of Ohm's Law which 
results when the left side of (2.55) vanishes. The modifications of (2.55) the 
generalized Ohm's Law which would ensue due to the influences of magnetic 
field, gravitation, pressure and viscosity would follow from (2.54) but are not 
required in the present discussion. 

One finds that the electrical charge density naturally decays (with or 
without oscillation) unless there is an ionic generation term which produces 
inhomogeneities in the differential equation for the charge density. These act 
as forcing terms and no matter how small they may be initially, they make 
possible the development of "exponentially increasing" charge and current 
densities. 

In section 3 a separation of variables argument is pursued to determine 
the relationships between the time-dependences of the charge density, current 
density and electric field without dropping the nonlinear term of the 
generalized Ohm's Law. It is found that there is a relative saturation of the 
convective current which allows the electric field to increase monotonically 
with time corresponding to a non-oscillatory displacement current. This 
continues until dielectric breakdown occurs which is accompanied by the 
development of an ionization channel with a greater conductivity than the pre- 
stroke medium. When this channel is developed a larger convective current can 
flow thereby relieving the high field which led to breakdown. After the 
channel has been developed the stochastic nature of the charge exchange process 
can lead to a succession of return and direct strokes. 

If one considers a crude model for field propagation along the ionization 
channel with the rather drastic assumptions that appreciable values of the 
fields are restricted to a cylindrical portion of the channel of radius r and 
length h and that all losses from this portion are radiative, one obtains for 
the power loss P 
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( 6 . 1 ) 


2 

P ■ S 2urh + S nr 
r z 

where S , S are the radial and axial components of the Poynting vector. If 
the power wire then taken to be produced by Joule heating one would have 

E J - 2(S /r) + (S /h) (6.2) 

z z r z 

For S w S and h >> r this would yield 

r z | 

r - (2 S /E J ) (6.3) 

r z z 

for the radius of the ionization channel. S could be obtained from the 

luminosity of the lightning flash and using r J - I/irr 2 (6.3) becomes 

z 

r - (2 S /irE n I) 1/3 (6.4) 

r d 

where E has been replaced by the breakdown field E 0 and I is the 
current. 

In the fourth section the statistical descriptions of the two charge 
species is given quantitatively in terms of the corresponding absolute 
densities of positive and negative charge which separate out in a thundercloud 
with the development of the electrical dipole moment of the cloud. In the 
fifth section physical arguments bearing on the charge separation processes are 
reviewed. 

It is believed by the author that considerable progress has been made in 
making the quantitative arguments of the first, third, and fourth sections with 
the work of the third section on discharge initiation being founded on the 
plasma-theoretic derivation of the generalized Ohm's Law in the second section. 
There is every indication that the concepts of section 5 can be supported by 
precise treatments of the Boltzmann transport equation in a subsequent 
publication. Finally a comment about triggering will be made. In addition to 
the field enhancement factors associated with excessive curvatures of portions 
of an aircraft one can expect triggering to occur because the aircraft 
themselves lower the impedance path for a lightning stroke. For metal-bodied 
aircraft this is because they decrease the free air path of the stroke and for 
composite-bodied aircraft because they raise the permittivity of the stroke 
locus. 
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